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ABSTRACT 


This thesis is centered around the construction and analysis of the principal 
arithmetic surface (3,5) over Q. By adjoining the two symbols i,j, where 7? = 3, 7? = 5, 
such that i7 = —j72, I can produce a quaternion algebra over Q. I use this quaternion 
algebra to find a discrete subgroup of SL2(R), which I identify with isometries of the 
hyperbolic plane. From this quaternion algebra, I produce a large list of matrices and 
apply them via Mobius transformations to the point (0,2), which is the center of my 
Dirichlet domain. This list of transformed points is my orbit sampling. The possible 
walls of the Dirichlet domain are perpendicular bisectors of the hyperbolic geodesic seg- 
ments between (0,2) and points from my orbit sampling. Once I produced the list of all 
perpendicular bisectors, I plot them. From this plot, I refine the collection to only the 
walls for my Dirichlet domain. I then analyze this surface by finding its hyperbolic area, 
and the matrices that correspond to gluing one side to the other. 

The production of this hyperbolic surface by adjoining (3,5) over Q is an “ex- 
ample problem”. This example problem produces a hyperbolic polygon that could be 
used to show others the process of finding a Dirichlet domain, or producing a hyperbolic 
polygon. One reason why one may want to find a hyperbolic surface is for the list of 
generators that correlates to the edge gluings. This list of generators can be used in other 
areas of mathematics and science. 

This project brought together many areas of math, including topology, hyper- 
bolic geometry, linear algebra, group theory, and noncommutative algebra. Additionally 
this project involved a significant programming component, in which I wrote and imple- 


mented code in Sage, a Python based computer algebra system. 
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Chapter 1 


Introduction 


In Larry Guth’s Metaphors in Systolic Geometry, he states “Arithmetic Hyper- 
bolic surfaces are among the strangest and most interesting examples in (Riemannian) 
geometry”. I decided to study hyperbolic geometry because I thought it was different. 
What made it different was that this was the first type of non-Euclidean geometry that I 
have experienced. I remember the first day learning about H?, the hyperbolic plane. My 
mind is blown when learning that in Euclidean geometry, we have lines but in hyperbolic 
geometry we have circles that are geodesics. When learning about surfaces, I found it 
was very hard for me to visualize the concept of bending and stretching of surfaces. This 
actually motivates me to build the surfaces I am working on by hand with clay. It gives 
me a deeper understanding to construct what covers a genus 2 or a genus 3 surface. I use 
straws as a “literal pole” (see Figure 2.14), to understand how the fundamental domain 


is being rotated. 


Figure 1.1: Clay figures-1 


Figure 1.2: Clay figures-2 


Once I had these constructions, it was the quantitative features of these surfaces 


that intrigue me the most. For example, two surfaces could look the same, but they could 


have different cone points (see Section 2). A hyperbolic surface can be understood to be 
a hyperbolic polygon with gluing maps. A special class of such polygons can be built 
from arithmetic objects called quaternion algebras. For my research, I create a hyperbolic 
polygon from the quaternion algebra (3,5) over Q. This algebra is formed by algebraically 
adjoining the square roots 3 and the 5 in a special way. More information on this 
quaternion algebra can be found in Chapter 4.2. 

The groups formed by this type of constructions are a special type of what are 
called Fuchsian groups. These groups have deep connections through out geometry and 
algebra. For example, in Small Generators of Cocompact Arithmetic Fuchsian Groups, 
by Michelle Chu and Han Li [Lil6], they state there is a special connection of the subclass 
formed from arithmetic Fuchsian groups to number theory and to automorphic forms. An 
automorphic form is a function from a topological group G to the complex numbers. The 
action for the discrete subgroup [TC G the topological group is invariant. Automorphic 
forms were first discovered as a generalization of trionometric and elliptic functions by 
the French mathematician Henri Poincare. For more information on automorphic forms, 
you can use the following resource [Deil2]. 

It is a nontrivial problem to determine a set of finite generators for Fuchsian 
groups. To determine the set of generators, you can create a special polygon in H? called 
a fundamental domain and the generators then correspond to certain side pairings of 
the polygon. The polygonal fundamental domain can be found in certain contexts using 
algorithms by Johansson, Voight and Macasieb, and others have written computer codes 
to implement this strategy. I do the same. I use Sage, a computer algebra system. 

To create this polygonal fundamental domain, I first find the quaternion algebra, 
and use the quaternion algebra to produce matrices (See Chapter 4.2). I use the matrices 
to transform the point (0,2), which is the center of my Dirichlet domain (see Chapter 
3.4). The list of points which I produce from the transformation of (0,2) is my orbit 
sampling. Using the points from my obit sampling, I find the perpendicular bisectors 
between (0,2) and each point in my orbit sampling. I use the perpendicular bisectors to 
find my Dirichlet domain (see Figure 5.14). 

Once I have the Dirichlet domain, to check that I have the correct walls, I use 
three different area formulas: Gauss-Bonnet Theorem 2.25, the Arithmetic Area Theorem 


4.14, and Riemann-Hurwitz Theorem 2.29, to confirm the correct area of my founded 


hyperbolic polygon. The Arithmetic Area Formula is actually being used in many papers 
like [Mac08]. The Riemann-Hurwitz Theorem is named after the German mathematicians 
Georg Friedrich Bernhard Riemann and Adolf Hurwitz. The Gauss-Bonnet Theorem 
(Theorem 2.25) is actually named after two mathematicians Carl Friedrich Gauss, who 
has developed a version of the theorem, but has never published it, and Pierre Ossian 
Bonnet who has published a special case of the theorem. 

In finding this hyperbolic polygon, I will then have produced a list of finite 
generators for my Fuchsian group. This list of finite generators, could be used in other 


areas such as spectral geometry, systollic geometry, and many other fields. 


Chapter 2 


Surfaces 


2.1 Classification of Surfaces 


I now begin with some background on the topology of surfaces. For a detailed 
point-set topology reference, see [Mun18]. A surface is a shape that locally looks like the 
plane. For example, the Earth is a surface and it is locally flat, and it is for this reason 


people thought it was flat for so long. 


Definition 2.1. A (topological) surface is a Hausdorff topological space with a countable 


basis, each point of which has a neighborhood homeomorphic to an open subset of R?. 


Definition 2.2. Let S be a compact surface. An orientation on S is a continuous 


choice of an “outside” to the surface and S is orientable if it admits an orientation. 


Definition 2.3. An orbisurface is a Hausdorff topological space with a countable basis 
{U;} such that to each U;, there exists an open set U; such that U; & U;/T; C R? and 
I; C Isom(R?) is a finite subgroup. If 1; is a rotation group, then the image of (0,0) in 


U; is a cone point. The order of 1; 1s called to degree of cone point, and is denoted e;. 


Definition 2.4. A good orbisurface is an orbisurface which is the quotient of a surface 


under the action of a discrete group of homeomorphisms. 


When working with the topology of surfaces, we can bend, fold, and stretch 


them as long as we keep the orientation of the surface. We can think of surfaces like clay. 


Definition 2.5. Let S be a compact orientable surface. The genus of S is the number 


of “donut holes”. (A more rigorous definition can be found in [Mun18, Chapter 12/) 


Definition 2.6. A punctured surface is a compact surface with p many “marked 
points”, which we consider having been deleted. We denote the punctured orientable 


surface with genus g and p punctures as Sg p. 


Definition 2.7. The signature of a good orbisurface with genus g, p punctures, and r 


many cone points of orders €1,...,€r, 18 (g; C1, -.-€r3 P)- 


Theorem 2.8. Good orientable orbisurfaces are completely characterized by 8 pieces of 


data: genus g, punctures p, and cone points (with orders) (e1...e,). [Mar12, Theorem 1.1] 


In other words, if two such orbisurfaces have the same signature, then they are 


homeomorphic orbifolds. 


2.2 Group Actions and Covering Spaces 


Definition 2.9. A discrete group is a group with discrete topology. 


Definition 2.10. If X is a topological space and I is a group acting discretely on X, 


then a fundamental domain for the action of I is a closed subset DC X such that 


1. Uyer(yD) = X 


2. Ify1 #2, then yD°172D° = 6. 


Definition 2.11. A covering space of a space X is a space x together with a map 
Dp: ano ¢ satisfying the following condition: There exists an open cover {U,} of X such 
that for each a, p~!(Uq) is a disjoint union of open sets in x: each of which is mapped 
by p homeomorphically onto U.. We do not require p~!(Uq) to be nonempty, so p need 


not to be surjective. 


Definition 2.12. A covering p: X — X is s-sheeted coverings if p is surjective and 
there exist an open cover of X, {Ua}, where p~!(Uq) is a disjoint union of s-many sets 


homeomorphic to Uq. 


The following is a well known result regarding orbifolds. It follows from the 


work of [Thu97, Chapter 13] and is useful in my research. 


Proposition 2.13. [fT is a finite symmetry group acting on an orbisurface, then this 
gives a gluing map on the fundamental domain. This newly glued surface is the quotient 


and we get a covering map with |T|-sheets. 


I start with building an understanding of the topology of surfaces. The topology 
of surfaces refers to the qualitative characteristics of a surface. When a surface covers 
another, this brings the opportunity for a surface to have cone points. Imagine taking 
a literal pole, see Figure 2.14, and intersecting this pole with a surface, a fundamental 
domain is rotated around this pole at a certain degree, covering the surface. The points 


of intersection between the pole and the surface are where the cone points are made. 
Example 2.14. Let’s look at an example of a genus 2 surface. 


Figure 2.1: A of genus 2 surface covering genus 1 surface with 2 cone points 
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Here we are starting with a genus 2 surface. By placing this pole in the middle 


of this surface, we have two points of intersection between the pole and the surface. That 


is the top entry point and the bottom exit point. We rotate the surface 180° around this 
aris. The right side (green) is a fundamental domain. Since our fundamental region is 
rotated by 180°, our two intersection points in the quotient will be cone points with 180° 
around them. So in this erample, a genus 2 surface covers a genus 1 surface that has two 


cone points that both have 180° around them. This is a 2-sheeted cover. 


Going forward, in these context, I will refer to this pole around which we are 
rotating the surface as the axis of rotation. Looking at the example above, let’s notice 
where I place the axis of rotation in the figure. It is directly in the middle between the 
two genus. You can place the axis in more than one place. The fundamental domain 
must be symmetric to the surrounding region. Let’s look at an example where I place 


this axis in two different places. 


Example 2.15. Here I look at two distinct axes of rotation for a genus 3 surface, and 


analyze the resulting quotient orbifolds. 


Figure 2.2: Genus 3 surface with two different placements of the axis of rotation 
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One surface has a fundamental domain that is being rotated by 180°. This 
creates four cone points that each have a degree of 2. In the visual, where the pole is 
placed directly in the middle, the fundamental domain is rotated by 120°. The two cone 
points both have a degree of 38. The way to distinguish between surfaces is by the signature. 
A genus 8 surfaces covers two genus one surfaces with the signatures of (1; 2, 2, 2, 2) 


and (1; 8, 3). There exist a 2-sheeted cover ¢, : (3;—) — (1;2,2,2,2) and a 3-sheeted 
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cover 9; (3; —) > (133, 3). 


The following result which is useful in my work is a well known theorem regarding 


surfaces. It follows from the work of [Mun18, Chapter 12]. 


Theorem 2.16. A genus g surface can be formed from a 4g-gon with appropriate edge 


gluings. 
Proposition 2.17. A genus 2n +1 surface, where n > 1, covers a genus 3 surface. 


Proof. Start with a genus 2n + 1 surface. Let’s first start by placing one of the holes in 
the center and then the others circled around it. Looking at the figure below, we can 
group together the adjacent pairs of holes around the center hole. The individual pair 


along with the center hole is the fundamental domain (see Figure 2.3). 


Figure 2.3: A genus 2n + 1 covers a genus 3 surface 
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By placing an axis of rotation through the center hole, you can rotate the fun- 


damental domain (22°. Upon identifying the boundaries of the fundamental domain, 


you get a genus 3 surface. 


Proposition 2.18. A genus n-+1 surface for alln > 1 covers a genus 2 surface. 


Proof. If we start with a genus, (n+ 1) surface, we can place one hole in the center, then 
the others around it. Notice by putting one hole in the center of our figure, we have n 
surrounding it. We can then group together each hole with the center hole. This group of 


paired up holes is actually the fundamental domain. Now by placing the axis of rotation 
360° 


through the center hole, we can then rotate the fundamental domain by ——. So a genus 


n+ 1 surface covers a genus 2 surface. 
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Figure 2.4: A genus n+ 1 surface covers a genus 2 surface 
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The geometry of surfaces is the “quantitative” features of shape. To learn about 
the geometry of surfaces, we dive into the hyperbolic plane, the hyperbolic metric, and 


hyperbolic surfaces. 
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2.3. Hyperbolic Geometry 


The upper half plane of R? is a model for the hyperbolic plane. Listed below 


are some general definitions that can be found in [Kat92]. 


Definition 2.19. The hyperbolic plane, H?, is modelled by the upper half-plane {(x, y) 


2 2 
€ R?|y > 0} endowed with the hyperbolic metric, ds? = oa 


Definition 2.20. Jf c: [a,b] > H? andt — (x(t), y(t)) is a curve, its length is €(c) = 


_ ph Ve)? +(y'()? 


a 


Definition 2.21. A geodesic is a curve that locally minimized length. 


Theorem 2.22. The geodesics in H? are lines and semicircles which all are perpendicular 


to the x-azis. [Buc11, Theorem 2.7] 


Definition 2.23. In the upper half-plane model of hyperbolic geometry, the area of a 


region R described in Cartesian coordinates, denoted as A(R) is given by 


A(R) = | if ately 


Definition 2.24. A hyperbolic polygon is a polygon in H? U OH?, where OH? is the 
x-axis, whose edges are geodesics. If at least one of it vertices lie on the boundary, we 


say it is an tdeal polygon. 


In order to find the area of hyperbolic surfaces, we use hyperbolic triangles. 
Hyperbolic triangles are like Euclidean triangles, instead of three points being connected 
by line segments, the three points are connected by geodesics. We do this by using the 


following Gauss-Bonnet theorem: 


Theorem 2.25 (Gauss-Bonnet). Jf a triangle is a hyperbolic triangle with angles a, 3,7, 
then the hyperbolic area of the hyperbolic triangle is 7 —a—B-—+y 


Proof. If a triangle has angles a, 6,7, there are two possible cases. Case 1: If a triangle 
has the angles a, 6, such that y is an ideal point. We can use transformations such as 
Mobius transformation discussed in Chapter 3.1 to send the ideal point to infinity and 


send circle c to the unit circle (see Figure 2.5). To find the area of the shaded region T 


(see Figure 2.5), then we will use integration. 


cos 8 
HS dxdy _ i Saye 
—cosa J V1—2x2 y? 
cos 8 1 oo 
=) =] dx 
— cosa 7] 1-2 
fls)- ==) 
= dx 
— cosa ee) 1 — 2? 


T T 

= 5 a5 - (ra) 
T T 
ge ae 
=T-a-—P 


Figure 2.5: Triangle with ideal point 
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Case 2: The intersection of three geodesics gives us a hyperbolic triangle. If the 
triangle has no ideal points, then it looks like 7, which is the area of they hyperbolic 
triangle we want to find. By creating another geodesic, the red geodesic in Figure 2.6, 
we create triangle 7), which is an ideal triangle that is proved in case 1. We then use the 
area of 75 to help us find the area of T,. To find the area of triangle 7), we subtract the 
area(T; U T2) — area(T>) shown in Figure 2.6. 


area(T; UT,) =7-—o-—(64+7)-6 
area(Th) = 7 —o-—0(m—a) 


area(T,) = area(T, UT2) — area(T2) =17-—a-—B-y 


Figure 2.6: Area of triangle with no ideal points 


From this theorem, we have the following consequence: 


Corollary 2.26. If p is a hyperbolic n-gon with interior angles a1, ...,Q%p, then area(p) = 


(n — 2)m — a1 — AQ —---— An 


Proof. To find the area of a hyperbolic polygon, we split it into hyperbolic triangles. We 
do this because we know how to find the area of a hyperbolic triangle (see Theorem 2.25). 


We use two different strategies to do so. Say we have a polygon with n number of sides. 


1. Using 1 vertex the n-gon, we can draw line segments to the nonadjacent vertices. 
Doing this splits the polygon into n—2 triangles. Recalling Theorem 2.25, we know 
the area of a triangle is t—a—(—v7, where a, (, and ¥, are the angles of the triangle. 
To find the area of this hyperbolic polygon, we add the areas of the triangles. 
area(t,) + area(t2) + area(tz) + area(ta) + area(ts) + area(tg) +... + area(tn_2). 


17 


Now the area is area(polygon) = 7 - N-O-—E+n-—P-a7z+n-—A-ap—Gt+r- 
H-I-B4+n-C-—J—K+nr—-—D-L-M-... Buta, = A+B4+C4+D+E+4+F, 
a3 =G+H,ag=I+J,a5 =K +L, 06 =O+P, Am_2) =M+N. So we can 


rewrite this as area(polygon) = m(n— 2) — ay — a2 — a3 — 4 — 5 — 1G — «.-» — An—2 


Figure 2.7: strategy 1 


2. Using 1 point in the interior of the polygon, and drawing line segments from the 
interior point to each vertex, we split the polygon into n triangles. Recalling The- 
orem 2.25, we know the area of a triangle is 7 — a— 6-7, where a, 6, and 7 
are the angles of the triangle. Now the area(polygon) = area(T,) + area(T2) + 
area(T3) + area(T4) + area(Ts) + area(Tg) + area(T7) +... + area(Tn) So we have 
area(polygon) = x—-S—I-J+n—-V—-—K—-L+n—-W-—-M-—N+n-O-A-B+ 
an-P-—C-D4+a-Q-E-F+427-R-G-H4... But, ay = A+N, a2 = BHC, 
ag =D+E,a,=F4+G,0a5=J+K,a5=M+L, a, = H+]. We can rewrite 


this as area(polygon) = m(n — 2) — a1 — a2 — a3 — a4 — A5 — AG — «.-. — A(n—2) 
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Figure 2.8: strategy 2 


Proposition 2.27. The area for a hyperbolic surface of genus g is A = 4r(g —1). 


Proof. Suppose we have a hyperbolic surface of genus g with no cone points. We take 
surface g, and we can map it to its associated polygon. For example, a torus is a genus 1 
surface, and it maps to a rectangle. Every genus will be a 4g-gon. A point has 360° or 27 
radians, when opening a surface, points are ’split” to ”unglue” a surface. When a point 


is ”split” these points will add up to 27. From Corollary 2.26, (n — 2)m — ay — ... — An 


where n is the number of sides. The polygon will have 4g number of sides, so we will 


replace n with 4g. So we have A = (4g — 2)m — 27, and by factoring and adding like 


terms, we can rewrite this as, A = 4a(g — 1). 


Proposition 2.28. The LCM of the orders of cone points divide the number of s-sheeted 


covers. 


Proof. By Definition 2.4 there exist a smooth and compact surface S with boundary and 
H < Sym(S) the symmetry group of S, is finite such that S > S/H & A. Let p bea 
preimage of a cone point. We have that stabyp = {h € H|h(p) = p}. The stabyp < H 
is a subgroup where the size of H is s, which is the number of sheets, and the size of 


stabzyp is the degree e; of the cone point associated with p. Similarly, we find cone points 
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with the orders of €1,...,e€n. By Lagrange’s Theorem (group theory) we know the order 
of the subgroup stabyp divides the order of H, therefore e;|s. Since e;|s for each i, 


LCM (e,--:en)|s. 


Understanding surfaces as polygons with edge gluings, I am now able to compute 


the area of hyperbolic surfaces. 


Theorem 2.29 (Riemann-Hurwitz Theorem). Jf S is a hyperbolic orbisurface with sig- 
nature (g;€1,..-,€n), and cuts open to a hyperbolic polygon P, then 


n 


area(S) := area(P) = 27 (2 24 Ss" a ‘ , 


cs 
i=l ¢ 


with n-many cone points e1,...,@n, let i denote indices which range between 1 and n, and 


g is the number of genus a surface has. 
This formula outputs the area of surface given its genus and cone points. 


Proof. Let’s first examine the area of a surface with 0 genus, and n cone points. We can 
open this surface to a polygon with 2n many edges. We first start by cutting this surface, 
and by opening it up we can find its polygon. Now using Corollary 2.26, we split this 


polygon into triangles. 


20 


Figure 2.9: Area of surface with n many cone points. 


Cn £4 


Using Corollary 2.26 we know that the area of a polygon with n cone points is 


(n — 2)m — a, —-++:— Qn. Now to find the area of surface with n many cone points and 
one boundary. Let’s say we have a surface with 0 genus and cone points ej, e2,e3 with 
one boundary. We can create a copy of our surface, and paste them together to have 1 
complete surface (see Figure 2.10). We know how to find the area of this surface now that 
it is whole. By dividing this area by 2 we compute the area of a surface with n many cone 
points and one boundary. Therefore, area = a qt ten)) which will simplify 


to area = 2a[(n —1) —(44+... ++). 


el en 
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Figure 2.10: Area of surface with n many cone points and one boundary 


es ft, 


Now consider a genus g surface with one boundary component. By Proposition 
2.27 the area of a genus g surface given is area(genus) = 4n(g — 1). Now using the 
same process above we have the area of a surface with genus g and one boundary is 


An(2g—1) 
2 


area(genus) = and this will simplify to area(genus) = 27(2g — 1) (see Figure 


2A iy 


Figure 2.11: Area of surface with g genus and one boundary. 


boundary 
Giving 


Therefore the area of a surface with g genus with n many cone points is Area = 


2n(2g —2+n— (2 ++ ~)) (see Figure 2.12). 
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Figure 2.12: Area of surface with g many genus and n many cone points. 


cut 
cunt 4 —> irenv 


Using this formula, I am able to write a Sage code that gives me my own list 
of surfaces that could possibly be covered by a genus 2 surface. Creating my own list 
using Sage and hand drawn constructions gives me motivation to then produce a list of 


surfaces that are covered by a genus 3 surface. I go over this in the next section. 


2.4 Genus 2 Surface Coverings 


I begin to study what is covered by a genus 2 surface. This leads me into 


studying Lemma 3.2 from [Mac08]: 


Lemma 2.30. /Mac08, Lemma 3.2] Let S be a hyperbolic orbisurface covered by a genus 
2 surface. Then S has one of the following signatures: (0; 2, 8, 7), (0; 2, 3, 8), (0; 2, 8, 
O), (08 2 28, 10), (02 2.2, 12); (0s 2, 2, 8)y (022, a, 6), (0 2 8) (Oe 2, ay 12), (02 2, 
Bi) (02 2.9.0), (0s 2, 0, 10), (05 2,6, 6), (0. 2, 8 BY. (0s 8, By A), (Op By oy BO) A Oe 8 
by 0),.. (0s 8; 8 8), (08 Bs, 2), (O08, 05.6), (OP ai 4 oy (02S, 2, 8), CE 2, 2 2 8), (0: 
Be Be ae (0) Be 2, By O)s (0) 2 2S a), (OS 2d gi (Ora, oy oo), on (0 2, 2, 2: 
2.2), (0: 2, 2, 2,2, 2, 2), (12 2), (1s 8), of (1r 2, 2). 


a 


This result claims to produce a list that contains all orbisurfaces which a genus 
2 surface covers. While examining [Mac08, Lemma 3.2], I produce constructions of or- 
bisurfaces which the genus 2 surface covers. In doing this, I have found one surface that is 
not on Macasieb’s list. However, it is unknown to me how to construct all of the surfaces 


that Macasieb claims are covered by a genus 2 surface. 
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In the following figures, I show the progression of surface coverings for a genus 
2 surface. Here is an example of how read the information in the figures below. Starting 
with a genus 2 surface, there are three ways to place the axis of rotation. The three ways 
to place the axis are numbered as 1, 2, and 3. Let’s take a look at 1. If I place the axis 
of rotation in the middle, then the left and right side of the genus 2 surface is symmetric. 
Notice above the axis of rotation there is the degree of rotation. Let’s look at the shaded 
region on the right side of this genus 2 surface, this is the fundamental domain. For 
this example, the fundamental domain is rotated 180°. On the right hand side of the 
genus 2 surface we started with you can see the quotient with its associated signature. 
Notice that this new surface (1;2,2) still has symmetry. For this specific surface, we can 
place the axis of rotation in two different places. These two different axis of rotation 
are indexed as l.a and 1.b. Let’s look at one way to find the covering for the indexed 
surface l.a. The image on 1.a, shows the orbisurface with cone points, as well as the 
rotation transformation of 180°. This rotation transformation produces the signature (0; 
2,2,2,2,2). The surface corresponding to this signature is still symmetric. I have found 4 
ways to place the axis of rotation on this surface, it is indexed as l.a.1, 1.a.2, 1.a.3, and 
1.a.4. I use the indexing this way to show you this was the surface found from the index 
l.a. 

Let’s look at 1.b. We start with the same signature as 1.a, (1;2,2), but I place the 
axis of rotation in a different location. 1.b produces the signature (1;2). Notice this surface 
still has symmetry, and it produces the signature (0;2,2,2,4), but this surface has already 
been covered above, so I write the indexing of where you can find the corresponding 
surface coverings for this surface. 

If the letters in the indexing have the same “place value” that means it is the 
same figure with the axis of rotation being placed through the surface in a different 
location. If that surface has symmetry, then it is continued, and I add the number 1 
following the letter to show that it is a continue to the previously used axis of rotation. 


The indexing is continued in this form. 


Proposition 2.31. For a genus 2 surface, I have constructed all of the following signa- 
tures: (0;2,8,8), (0;2,3,9), (0;2,8,12), (2,3,18), (0;2,4,5), (0,2,4,6), (0;2,4,8), (0;2,5,10), 
(0:2,6;6), (0;2,8;8), (0;3,3;4), (0:3, 454), (0:3, 6,6), (0:4,454), (0:2,2,2,8); (0;2,2,2,4),. (032, 
29.3), (022.44), (0:3;3)3,3)) (O2,2,2,22); (O22 2,2 2.2), (132), (2,2). 
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Proof. Proof by construction. 
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iv 
(0, 2.4, b) 
40° 
. (9, 3.9.4) 
; (0, 1.3.19) 
120° 110° 180° (0 


While constructing these signatures, I have found the signature (0; 2,3, 18) that 
is not listed in Macasieb’s Lemma 3.2. There are signatures that I can not construct 


on Macasieb’s list, so I further my research to see if either I can rule out any signatures 
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or add to my list of constructions. In doing this I also want to confirm my findings 
of the signature (0;2,3,18). To do this, I use Sage coding to help in my research and 


calculations. 


2.4.1 Sage Coding - Producing Signatures of Surfaces that Cover a 


Genus 2 Surface 


Sage is a Python-based programming language. I use Sage coding throughout 
my research. Using the Riemann-Hurwitz Theorem (Theorem 2.29) and Proposition 2.28, 
I create a Sage code that produces possible signatures for a genus two surface. Since we 
start with a genus 2 surface, the surfaces that cover a genus 2 are a genus 1 and genus 
0. This constrains possible number of cone points and degrees. From this, I produce a 
finite list of signatures that may be realizable. In this code, we find a range of values for 
cone points. From Proposition 3.23, we know the area of the quotient orbisurface divides 
4r(g — 1). For a genus 2 surface, the quotient must divide 47. To find the maximum 
number of cone points reasoning about the fact that when covering a surface, the quotient 
often has higher cone points but a smaller amount of cone points. The greatest number 
of possible cone points comes when the surface has cone points of smallest degree which is 
order 2. When covering a surface, there has to be symmetry to rotate the surface around 
the axis of rotation. The smallest order is 2 which is 180°. We could use the function 
4n = (interger > 1) -2a(—2+1r(1/2)), where r is the number of cone points. In solving 
we get there are 6 possible cone points for order 2. This is the process I go through to 
find possible signatures. 
for a in range(2,3): 

for b in range(a,7): 

M=[2,a,b] 

L=1cm(M) 

temp=0 

for m in M: 

temp += (m-1)/m 


if temp<=(2+2/L): 
print((2,a,b)) 


for a in range(2,3): 
for b in range(a,7): 
for c in range(b,100): 
M=[2,a,b,c] 
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L=1cm(M) 
temp=0 
for m in M: 
temp += (m-1)/m 
if temp<=(2+2/L): 
print((2,a,b,c)) 


for a in range(2,10): 
for b in range(a,10): 
for c in range(b,100): 
for d in range(c,100): 
for e in range(d,100): 

M=[a,b,c,d,e] 

L=l1cm(M) 

temp=-2 

for m in M: 
temp += (m-1)/m 

if temp<=(4/L) and temp>0: 
print ((a,b,c,d,e)) 


These codes are examples of how possible signatures can be found. Shown above is a 
genus zero surface with 3, 4, or 5 cone points. Macasieb’s paper [Mac08, Lemma 3.2] has 
a smaller list. It is an open question as to why she excludes certain signatures. With this 
code, I notice that all of the signatures listed on Macasieb’s paper are on mine except the 
signature (0; 2,5,6). Iam able to exclude this signature from my list using the Riemann- 
Hurwitz Theorem (Theorem 2.29). Macasieb does not list the signature (0; 2,3,18), and 


I construct it in Proposition 2.31, and also produce it using Sage. 


2.5 Genus 3 Surface Coverings 


Finding the constructions of quotients of a genus 2 surface motivates me to find 
the constructions of a genus 3 surface. In the figures below, I show the progression of 
surface coverings for a genus 3 surface. Each image shows the orbisurface with cone points, 
as well as the degree of the rotation transformation. Each surface on the right hand side 
of the figures shows an arrow to the right, that arrow points to the quotient as well as the 
associated signature. The indexing is the same as Proposition 2.31 constructions. Notice 
that a genus 2 surface covers a genus 3 surface, and the constructions for the genus 2 


surface is shown in Proposition 2.31. I list them here, but I do not reconstruct them. 
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Proposition 2.32. For a genus 3 surface, I construct all of the following signatures: 
(0;2,8,8), (0;2,8,9), (0;2,3,12), (2,8,18), (0;2,4,5), (0,2,4,6), (0;2,4,8), (0;2,4,12), (0;2, 
5,0) (0;2,9,10); (0:2,6;6) 5 (O22,6,9) < (O32, 7,12): (072,6,8), (02,912.12), (0:3, 3,4), (03,2, 
4);(0;8,0,0), (04,44): (0:4,4,0), (0;2,2,2,8), (0;2,2,22)) (02,226); (02, 2,8:3), (02:2, 
4,4), (052,260); (03,3008), (0:2:2,9,9,8), (02,22 22), (022,242), (0:2,2,2,2 2:2), 
(02 ZOD 2 OLD) (OG heh TiVo) (igo (lee 2.2) and (250). 


Proof. Proof by construction. 
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The above constructions show what a genus 3 surface covers. Notice that a 
genus 2 surface covers a genus 3 surface, see Proposition 2.31 for genus 2 constructions. I 


will now use Sage coding to see if there could be any signatures that I could not construct. 
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2.5.1 Sage Coding - Producing Signatures of Surfaces that Cover a 


Genus 3 Surface 


Using Sage, I created code that will produce a list of signatures which can be 
covered by a genus 3 surface. This code is created using the same reasoning as Section 
2.4.1, but a genus 3 surface can be covered by a genus 2, genus 1, and a genus 0. 


for a in range(2,10): 
for b in range(a,10): 
for c in range(b,100): 
M=[a,b,c] 
L=1cm(M) 
temp=-2 
for m in M: 
temp += (m-1)/m 
if temp<=(4/L) and temp>0: 
print ((a,b,c)) 


for a in range(2,10): 
for b in range(a,10): 
for c in range(b,100): 
for d in range(c,100): 

M=[a,b,c,d] 

L=l1cm(M) 

temp=-2 

for m in M: 
temp += (m-1)/m 

if temp<=(4/L) and temp>0: 
print ((a,b,c,d)) 


for a in range(2,10): 
for b in range(a,10): 
for c in range(b,100): 
for d in range(c,100): 
for e in range(d,100): 

M=[a,b,c,d,e] 

L=l1cm(M) 

temp=-2 

for m in M: 
temp += (m-1)/m 

if temp<=(4/L) and temp>0: 
print ((a,b,c,d,e)) 


This Sage code produces the possible signatures of a genus 0 surface that has 3, 4, or 5 


cone points. 
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Proposition 2.33. Using the Riemman-Hurwitz Theorem, for a genus 3 surface I produce 
the following list of possible surface coverings: (0; 2, 3, 7), (0; 2, 3, 8), (0; 2, 3, 9), (0; 
Do: TO). (ORE So Io) oe Be te (Oe Bo toe CON eo TO) LO ex 8,4)? Ms eo, 
D0) (On 224s DS). A Os-2 hy 0), (Os By de. TV), (Or 2, he 8) (02 Og 10) (Os 2s 1B) Als D 
4,16); (0 2. de 20): (0 2: 5), (0; 2), (6), (02 25; 10), (022; 6,6), (0p 2; 69). 0; 
DOs 12 Os Bo De HOOD hee TAN TOFD. BoB). £08 EO). POR 8. SO) CO oe BoB) 
(ONS 8ST), (OP38, 8) D108. 8: SID) (08 By 8. 15). (OF By A, By (ON By de CO) (Oe 8s eg, 
12), (0; 3, 5, 5), (0; 3, 6, 6), (0; 3, 9, 9), (0; 4, 4, 4), (0; 4, 4, 6), (0; 4, 4, 8), (0; 4, 8, 
SB (OPS, 8,3 )s/ (02 6.66), (08T, TD (Oe ey 2: BOP, oO d) (O82 Oy 2, 8) 00% 
22 2D iO) ih Ot OO De OP, BP IO PO ee DO da (Ol ee, Opi On eu ee eds 
(0? DDO LC) (OSS, Bs OED Be Be Od, COE Dads ches Ad COPD, Bay Odie (Op Dy cley Aes 
4): (O72. See) Oe SS SS) (OR 2, Bo 2a), (O02 BOB 5 dg COS ae Soe a, 
Ds ( 0,222 Layee) 


I produce this list using the Riemann-Hurwitz Theorem and Sage, but there 
maybe signatures that cannot be realized. The signatures that are constructed in Propo- 
sition 2.32 are proven to be actual signatures for surfaces that are covered by a genus 3 
surface. It is an open question as to how to certify that the remaining ones on Proposition 


2.33 that I do not list in Proposition 2.32. 


Proposition 2.34. The list of genus 3 surface coverings found in Proposition 2.893, 
but can not be constructed in Proposition 2.82 is: (0; 2,3,7), (0;2,3,10), (0;2,3,14), 
(O:2,3,15), (0;2,3,24), (0;2,8,80), (0;2,4, 7); (0:24; 10); (02,4, 16), (0;2,4,20);, (0;2;5,6), 
(0;2;6,9), (02,612), (0:2,%, 7); (0;8,8,0),. (078,850); (053,38, 1)». (038;3,8) (078,512), (053; 
8,15), (0;3,4,6), (0;3,4,12), (0;3,5,5), (0;3,9,9), (0;4,4,8), (0;4,8,8), (0;5,5,5), (0;6,6,6), 
(O77), (02,2,2,9),. (0;2,2,2,8),. (0;2,2,2,10), (0;2,2,8,6),. (0;2:2,6,6),. (0;2;,8,3,3), (032, 
3,3,6), (0;2,4,4,4), (0;4,4,4,4), (0;3,3,3,3,3). 
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Chapter 3 


Hyperbolic Transformations 


3.1 Mobius Transformations 


A question I ask myself is: what is the difference between a surface and a 
hyperbolic surface? Surfaces are built from polygons along with their associated gluing 
maps. A hyperbolic surface is built from hyperbolic polygons along with their associated 
gluing maps. A hyperbolic polygon has edges that are geodesics. The gluing maps for 
hyperbolic polygons are isometries of H?. Listed below are general definitions that can 


be found in [Hat01]. 


Definition 3.1. An isometry of H? is a transformation T : H? — H?, that perserves 


the hyperbolic metric. 
In my research, I use Mobius transformations of H?, denoted as Méb(H?). 


Definition 3.2. A Mébius transformation is a transformation T : H? + H?, z > 


az-+b 
cz+d’ 


for any real numbers a,b,c,d such that ad — bc £ 0. 


Theorem 3.3. /Leh14, Proposition 5] Mobius transformations preserve angles. 
Theorem 3.4. /Kat92, Theorem 1.3.1] Méb(H?) © PSL(R) 


Theorem 3.5. /Bon55, Proposition 2.18] Mobius transformations send circles and lines 


to circles and lines. 


Let Isom(H?) be the set of all isometries, and I[som* (HI) be the subset of all 


orientation preserving isometries. 
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Theorem 3.6. /Bon55, Theorem 2.11] If T €Méb(H?) then T € Isom(H?). Further- 
more, Isom* (H?) ¥M6b(R?). 


An example of an isometry is the identity map, translations, rotations, and 
reflections. In my research, I use Mobius transformations for the edge gluings of my 
hyperbolic polygon. For convince, I work in SZ2(R) and I keep in mind, A“ =” — A of 
Méb(H?). 


Example 3.7. Mobius transformations are 


1 6 
1. 23> 2z+5 given by matrix , which is a parabolic transformation; 
0 1 


0 
2.24 (4) given by matrix , which is an elliptic transformation; 
1 


a 0 
3. 2 a®z given by matrix A this is a hyperbolic transformation. 
Oa 


Let A € SL2(R). 
Definition 3.8. The discriminate of A is disc(A) = b? — 4ac. 
Definition 3.9. The trace of A is tr(A) =a+4d. 


Proposition 3.10. /Kat92, Section 2] There are three fundamentally distinct types of 
Mobius transformations. Let A € SL2(R) and let T,4 denote the associated Mobius trans- 


formation. 
1. A is hyperbolic: 
e disc(A) >0 
e |tr(A)| > 2 
e A is diagonalizable 


e 74 has precisely two distinct fixed points in H?, both of which lie on the bound- 


ary. 


2. A is parabolic: 
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e disc(A) =0 


\tr(A)| < 2 
e A is triagonalizable 


T, has precisely one fixed point in H?, and it lies on the boundary. 


3. A is elliptic: 


e disc(A) <0 
\tr(A)| <2 


e A is rotationalizable 


T has precisely one fixed point in H?, and it lies in H? 


3.2 Circle Inversions 


For my research, I use M6bius transformations to find the perpendicular bisector 
between two points. In this process, I use two types of inversions, inverting a circle over 
another circle, and inverting a inscribed circle over the circle it is inscribed in. Let’s first 


start with what is an inversion. Below I list general definitions from [Nee97]. 


Definition 3.11. A circular inversion is a transformation where point A in the Carte- 


sian plane is transformed based on a circle a with radius k and center O such that 
OA: OA =k? 


Figure 3.1: circle inversion 


A' 


45 


Proposition 3.12. /Nee97, Theorem 3.2 (5)] If inversion in a circle is centered at q 


maps two points a and b to a’ and b’, then the triangles aqb and b'qa’ are similar. 


Figure 3.2: Similar triangles under circle inversion 


a! K 


Theorem 3.13. /Nee97, Theorem 3.2 (7)/ If a line L does not pass through the centre q 


of K, then inversion in K maps L to a circle that passes through q. 


Figure 3.3: Circle inversion maps circles to lines 


a’ L 


Theorem 3.14. /Nee97, Theorem 3.2 (9)/ If a circle C does not pass through the centre 
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q of K, then inversion in K maps C to another circle not passing through q. 


Figure 3.4: Circle inversion maps circle to circles 


Theorem 3.15. /Nee97, Theorem 3.2 (10)] Under inversion in K, every circle orthogonal 
to K is mapped to itself. 


Figure 3.5: Circle inversion circle maps to itself 


AT 


Theorem 3.16. /Kat92, Theorem 1.3.1] Circle inversion centered on the x-axis are 


orientati- on-reversing hyperbolic isometries. 


3.3. Discrete Subgroups 


We want to find the Dirichlet domain (see Section 3.4) because this give us my 
end goal of finding a hyperbolic polygon to build a hyperbolic surface. In order to find a 
Dirichlet domain, we use discrete subgroups of SL2R, which we are thinking of as gluing 
isometries of H?. In topology, a discrete subset of a space is one which you can separate 


all points with disjoint open balls. Below I list general definitions from [Kat92]. 


Definition 3.17. [fl C SL2R is a subgroup which is discrete relative to the subspace 


topology, it is a discrete subgroup. 


Definition 3.18. [fT is a group acting on a set X, and x € X, then the orbit of x is 
the set x = {g-a|g ET}, the full set of objects x is sent to under the action of Tx. 


Definition 3.19. The action of a subgroup, 1 C SL2(R) is properly discontinuous, 
if for any x € H?, and any closed and bounded K C H?, K AT =z is finite. 


Definition 3.20. A Fuchsian group is a discrete subgroup of PSL2(R). 


Proposition 3.21. [Kat92, Corollary 8.7] 1 is discrete if and only if its action on H? is 


properly discontinuous. 


Proposition 3.22. /Bucl1, Theorem 3.2] Fuchsian groups act properly discontinuously 
on the hyperbolic plane H?. 


Proposition 3.23. /Kat92, Section 3.6] Let T be a Fuchsian group, and F' be a funda- 
mental domain, then the quotient H?/T & F/T and the area(H?/T) := area(F). 


Discrete subgroups will help us to find the associated Dirichlet domain. 


3.4 Dirichlet Domains 


In order to find a hyperbolic polygon, I construct a Dirichlet domain. Dirichlet 
domains are an explicit process for producing fundamental domains for the action of 
discrete subgroups on H?. The following result, which is useful in my work, is a well 


known theorem on Dirichlet domains [Kat92, Section 3.2]. 
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Theorem 3.24. The process of finding a Dirichlet domain is: 
Step 1: Fix a point xo € H? to be the center of the domain; 


Step 2: Compute the orbit Tx; 


Step 8: For each x1 € Vx, find the perpendicular bisector between xox1; and 
Step 4: The domain, D, is the intersection of all half-planes containing xo. 


For Fuchsian groups, a Dirichlet domain only has finite sides. Therefore, I do not 
need the full orbit, just a well-chosen finite subset. This is fortunate, as it is impossible 


to construct the infinite orbit, but I can construct a well-chosen subset, which I will call 


my orbit sampling. 


Theorem 3.25. /Kat92, Theorem 11.8] Let {T;} be the subset of [ consisting of those 
elements which pair the sides of some fixed Dirichlet region F. Then {T;} is a set of 


generators for T. 
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Chapter 4 


Arithmetic Constructions 


4.1 Quadratic Extensions 


For my research, I produce a Dirichlet domain. To do this, we use a discrete sub- 
group. A discrete subgroup tessellates the plane with the Dirichlet domain, i.e. translates 
of out domain will not overlap. In order to get a discrete subgroup, we use quaternion 
algebras. To help in my understanding of quaternion algebras I have to gain an under- 
standing of quadratic extensions. Listed below are general definitions that can be found 
in [Voi97]. 


Definition 4.1. A quadratic extension of Q is the algebra A = Q|a] where a2 =n, n 


is square-free integer. 


Quadractic field extensions have embeddings into R and C. In order to know 
where the quadratic field extension is embedded, we look at the mapping. Does the poly- 
nomial have distinct roots, complex conjugate roots. For example, let Q{a] = Q[z]/(x? — 
2), where a + V2 and a > —V2. Q|a] has two real roots, both of which are embedded in 


R. Since we are working with quadratic field extensions, then we have two embeddings. 


Proposition 4.2. An element in Q|a] is of the form « + ya where x,y € Q. Q(a) x 
x ny 


C Mat2x2(Q) 
y 2 


On 
Proof. Let x,y € Q, where + ya € Qa] is an element. Let A = and t= 
1 O 
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1 0 On x ny 
So <a Sac +y = x,y € QP C Matox2(Q). 


1 O y x 


Theorem 4.3. If n is square free and x,y € Z such that they are primitive solution to 


2 


Pell’s equation, i.e. x? — ny? =1, then: 


xn 
1. The matrix corresponding to Pell’s equation d is a hyperbolic matriz. 


y 2 


2. The fixed points under a Mobius transformation are —\/n,/n 


zon 

Proof. Proof of 1: Let n be square free and x,y € Z. The trace of 4 is 2x. This 
y « 

matrix corresponds to Pell’s equation, where x = 1 is not a possible solution because if 


x = 1, then we would have —ny? = 0 and n must be a square free number. So x > 2, and 


xn 
tr Y > 2. By Proposition 3.10, this is a hyperbolic matrix. 
y «& 
Proof of 2: Let n square free and x,y € Z then 2? — ny? = 1. We have 
2 2 x ny . . oe . ° 
xz =1-—ny*. Then -r =r. Then it can be written as a Mobius transformation 
y «& 
(See Definition 3.2). So we have “*"" — r, Now we will solve for the roots, so xr + ny = 


yr +x 
yr? + ar. Now we will add like terms and have yr? — ny = 0. Now yr? = ny. So r = /n 


and r= —,/n. 


Definition 4.4. The ring of integers of a Q|a] is the subset R of Qla] consisting of 


elements that satisfy monic polynomials over Z. 


4.2 Quaternion Algebras 


Listed below are general definitions that can be found in [Voi97]. 


Definition 4.5. A quaternion algebra over Q is the algebra A = Qh, j,k], where 
a,b€ Qi? =a, j? =b, and ij = —ji. This is denoted with the symbol A = (a,b) 


For this project, restrict to the case a,b > 0. Quaternion algebras are like 
quadratic field extensions, it is the adjoining of two symbols, but they anti-commute. 
The maximal orders that come from quaternion algebras give us the discrete subgroup 


needed to find the Dirichlet domain. 
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Definition 4.6. The standard embedding I will use is p: At> SLo(Q|Va]) such that 
Je t Ja 01] 1|0 b 


1 b 
pilwh ,Jjr ,tbH , and k , and 
Oral 1 0 0 —-/a 0 -Va} }1 0 


extend linearly. 
Theorem 4.7. /Kat92, Theorem 5.2.7] If O C A is an order, and p: A > Mat2z2(R), 
then p(O) is discrete. 


We use the maximal order of © to then find the norm 1 elements, which will 


help us to find the quaternion algebras. Below is a list of definitions from [Voi97]. 


Definition 4.8. An order O C A is maximal if it is not properly contained in another 


order. 


Theorem 4.9. /Voi97, Proposition 2.3.1] A quaternion algebra (a,b) is either a division 


algebra or isomorphic to the matrix algebra of 2 x 2 matrices over R. 


Definition 4.10. [ft+ai+yj + zk € A, then the norm is: 


Nrm(t + 2i + yj + zk) = t? — ax? — by? + abz?. 


Theorem 4.11. /Kat92, Equation 5.2.11] If we have the embedding p: A + Mat212(Q|V/a)), 
anda € A, then Nrm(a) = detp(a) 


Theorem 4.12. /Kat92, Theorem 5.2.7] If A is a quaternion algebra (a,b) over Q, 
and O C A is a maximal order, and OW Cc O is the set of norm one elements, then 


T := p(OM) C SLo(R) is a discrete subgroup. 
Definition 4.13. A group of this form is a principal arithmetic Fuchsian subgroup. 


Theorem 4.14. /Bor81, Arithmetic Area Formula Theorem 7.8] If T is as in Theorem 
4.12, then the area of H?/T is: 


T 
area(H?/T) = Myeaalp 1), 


where AA is the set of primes for which A ramifies. 


Theorem 4.14 is very technical and powerful. The proof of this theorem is 
beyond the scope of the thesis. The Riemann-Hurwitz Theorem (Theorem 2.29) can be 
constructed by finding the area of a surface by breaking it into sections (See Figure 2.12). 
Unlike the Riemann-Hurwitz Theorem, the Arithmetic Area Formula is an arithmetic 
formula that does not require the qualitative aspects of a surface to find its area. These 


three area functions are from completely different perspectives. 
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Chapter 5 


The Research Process 


5.1 The Quaternion Algebra 


The goal of my research is to produce a hyperbolic surface associated with the 
maximal order in the quaternion division algebra (3,5) over Q. To do this, I find the 
Dirichlet domain. I center my Dirichlet domain at the point (0,2). I find the quater- 
nion algebra, which helps me to find matrices that transform the point (0,2). A list of 
sufficiently many transformed points is my orbit sampling. In order to find the orbit 
sampling, I first find the quaternion algebra by adjoining (3,5). I use Sage to help me 
find the maximal order, and the norm 1 elements. This produces a list of matrices. These 
matrices transform the point (0, 2). 


A.<i,j,k> = QuaternionAlgebra(3,5) 
O = A.maximal_order () 


#brute force ( Not Optimal) 
norm_one_elements = [] 
r = 10 
for a in range(-r,r): 
for b in range(-r,r): 
for c in range(-r,r): 
for d in range(-r,r): 
Oo = axwtb*x+c*yt+d*z 
if o.reduced_norm()==1: 
norm_one_elements. append (o) 


From the norm 1 elements, we produce matrices. 
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#turns norm one elements to a matrix 
def turn_order_to_matrix(o): 

R.<x> = PolynomialRing (QQ) 

K.<d> = NumberField(x*2-3) 

Id= matrix(K, 2,2,[1,0,0,1]) 

I= matrix(K, 2,2,[d,0,0,-d]) 

J= matrix(K, 2,2,[0,5,1,0]) 

K= Ix*J 

Matrix_Gens = [Id, I, J, kK] 

M = sum([list(o) [In] *Matrix_Gens[n] for n in range(4)]) 

return M 


def matrix_rounded(U): 
U_rounded = Matrix(2,2,[sigma(u) for u in U.list()]) 
return U_rounded 


Listed above, I have exact and rounded matrices. The exact matrices include the V3, V5. 
The rounded matrices, rounds the V3, V5, and matrix with 4 rounded decimal entries. 
Once we have the matrices produced from the function above, we use these matrices to 


create Mobius transformations. 


#take the founded matrices, and using the point (0,2) 
or complex number 21 to find where matrix takes 21. 
def mobius(B,PO0=2*sqrt(-1)): 
N1 = (B[0] [0] *P0+B[0] [1])/(BL1] [0] *P0+B[1] [1]) 
return N1 


5.2 Orbit 


Using the Mobius transformations found from the quaternion algebra we can 


now create an orbit sampling by applying these transformations to the point (0, 2) 


orbit = [J 
for o in norm_one_elements: 
B = turn_order_to_matrix(o) 
C = matrix_rounded(B) 
Pi = mobius(C) 
orbit .append(CC(P1) ) 
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Figure 5.1: Orbit sampling 


5.3 Finding the Dirichlet Domain 


To find the Dirichlet domain, I have to find the perpendicular bisectors of the 
geodesics that go through the point (0,2) and the points in my orbit sampling. 


Proposition 5.1. Given 3 noncolinear points, I can construct a circle through them. 


Proof. Let A = (x0, yo), B = (x1,y1), and A’ = (x%0,—yo) be 3 noncolinear points, by 


drawing line segments between them we can construct a triangle. 
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Figure 5.2: 3 points create a triangle 


B 
(x, ¥;) 


A 
(Xp; Yo) 


A! 
(Xp; Yo) 


The intersection of the Euclidean perpendicular bisectors for each edge is the 
circumcenter. To algebraically find the circumcenter, we create two of the euclidean 
perpendicular bisector lines and solve for the intersection. Let’s find the euclidean per- 
pendicular bisector between segment AB, let’s call the line that goes through points AB 
f(x), and let’s call the euclidean perpendicular bisector of AB g(x). To find g(x), we find 
the average the midpoint between points A, B. To do this, we take the average between 
the x and y values for points A,B, so we have a = 1 tto and b = iF yo | g(x) has a 
slope that is the opposite reciprocal of g(x). So the slope of f(x) is M = 42—“. Now 


©1—20 
the slope of g(x) is Mgz) = 57. 80 f(x) = m(a—a) +b and g(x) = (x — a) +b. Now 
we need another line to intersect g(x) with to solve for the circumcenter. Let’s find the 
Euclidean perpendicular bisector for the line that goes through the points AA’. Since 
line AA’ is a vertical line, the Euclidean midpoint is just the average of the y values, but 
since A’ = (xo, —yo), the midpoint is the x — axis. The Euclidean perpendicular bisector 


is y = 0. We set it equal to g(x) and solve for the x value of the intersection, which I call 


x2. In doing this, we have x2 = mb +a. 
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Figure 5.3: Production of the circumcenter 


(xy, yy) 


(Xp; Yo) 


x-axis (x, 0) 


To find the radius, we just find the distance from (x2,0) and point A. So 
r= / (x0 — £2)? + (yo)?. 


Figure 5.4: Circle through 3 noncolinear points 


Proposition 5.2. Given two points A,B € H?, we can construct a circle that goes 


through these two points. 


Proof. Let A,B € H? such that A = (x9, yo) and B = (x1,y1). There are four possible 
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cases we encounter in the creation of the circle, c, that goes through points A and B. 
Case 1: xp = x1 and yo = y1. This is the case where the points are the same. So there is 
no circle that goes through the two points. Case 2: 29 = 21. This is the case where the 
points are vertical. The geodesic that goes through these points is a vertical line that is 


centered at (%9,0) with a radius of oo. 


Figure 5.5: Line through points AB 


Case 3: yo = yi and yo = O. This is the case where the points are on the 
boundary. When this is the case, to find the center of the circle between points A, B we 
take the average between the x-values, I call this z-value a. The radius is the distance 


from the center of the circle to one of x-values. So the center of circle c is found from 


|(a — x0)|. 
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Figure 5.6: Circle AB on boundary 


boundary 


Case 4: This is the case where point A = (29, yo) and B = (21, y; are two distinct 
points. If two points are not equal, then by reflecting one of the two points over the x-axis, 


we can construct a unique circle through these three noncolinear points. By Proposition 


5.1 we create a circle centered at (%2,0) with a radius of r = \/(xo — £2)? + (yo)? (see 


Figure 5.4 for visual). 


This is shown in the code below. This code is classical and well known. 


def hyperbolic_geodesic_circle(P0,P1): 
#given two points, this function finds the geodesic 


x_0 = PO[O] 
y_O = PO[1] 
x_1 = P1i[0] 
y_1 = Pi[1] 


if y_0==y_1 and x_0==x_1: 
print("same point") 
#single equals is an assignment 
return( (0,0) ,0) 

elif x_0==x_1: 
print("same x’s but different y’s") 
r= "infinity" 
return ((x_0,0),r) 

elif y_0O==y_1 and y_0==0: 
print("both points on boundary") 
a = (x_1+x_0)/2 
center=(a,0) 
radius=abs (a-x_0) 
return (center,radius) 
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else: 
print("this is the generic case") 
a=(x_1+x_0)/2 
b=(y_1+y_0)/2 
m = (y_1-y_0)/(x_1-x_0) 
#f(x) = m*(x-a) + b #equation of line through starting points 
#e(x) = -1/m*(x-a) + b # 
x_2=m*bta 
r = sqrt ((x_0-x_2)72+(y_0)~2) 
return ((x_2,0),r) 


Figure 5.7 Below shows an example of the geodesic between the point (0,2), and (3,3). 


Figure 5.7: segment AB 


We will find the perpendicular bisector of segment AB. I create a function that 
uses a circle inversion to take this geodesic curve to a vertical line. To do this, I start by 
creating a circle, d, that goes through point A and has a center that is the left endpoint 
of circle C’. See the Figure 5.8 below. 
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Figure 5.8: Circle d 


To create circle d, we first find the left endpoints for circle c. To do this, I create 
a function that inputs the center and radius of the circle, and outputs the endpoints. Since 
a geodesic circle is perpendicular to the z-axis, our endpoints will be on the hyperbolic 


boundary. 


def circle_endpoints(P0,r): 


xO = PO[O] 
yO = PO[1] 
LO = x0 - r 
RO = x0 +r 


return ((LO,0),(RO,0),r) 


I use the left endpoint of circle c. 


def radius_of_circle_d(P0O,P1): #radius of circle d 


xO = PO[0] 
yO = PO[1] 
xi = P1i[o] 
yi = P1i[1] 


ri = sqrt((x0 - x1)°*2 + (yO - y1)*2) 
return (P0O,r1) 


The purpose in the creation of circle d is when I compute a circle inversion 


transformation of circle c being reflected over circle d, this moves circle c to the y-axis. 
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Proposition 5.3. If there is a hyperbolic circle inversion fixing A that takes circle c to 
the vertical line € going through A, then point B on circle c goes to point D, where D is 


the intersection of € with the vertical axis. 


Proof. By Proposition 5.2, I compute the center and radius of c which I denote center = 
(P[0], P[{1]), radius = r. To find where point B has moved, we create a line that goes 
through the center of circle d, and point B. The intersection of this line along with the 
y-axis gives us the location of the image of B under this transformation. In Figure 5.9 
below, this is point D. Notice that segment AB has moved to segment AD. Using the 
left endpoint of circle c which I denote as L where L = P[0] — r, we create a circle, d, 
that goes through point A. The inversion of circle c over circle d will move circle c to 
a vertical line that goes through points AA’ (See Theorem 3.5). By Definition 3.11 to 
find where B has moved, we can draw a line from L to point B. The intersection of line 
through ZB which I will call f(a), and the line through AA’ which I will call g(x), is 


the transformed B, called point D in Figure 5.9. The slope for f(a) ism = ant and 


using point B, we use the point slope formula to create a line. So f(x) = m(x—21)+ 1, 
now g(z) is the line through AA’, but lets recall from Proposition 5.2 that A and A’ both 


have the x-value of xp. so g(x) = 20, and the intersection is y2 = m(ao — 21) + yi. We 


have D = (20, y2). 


Figure 5.9: Segment AD 


I create a Sage function, that computes Proposition 5.3, the circle inversion that 


transforms point B to point D shown in Figure 5.9. 
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def circle_inversioni(P0O,P1): 
(P,r) = hyperbolic_geodesic_circle(P0O,P1) 


x_0 = PO[O] 

y_O = PO[1] 

x_1 = P1[0] 

y_1 = P1[1] 

if r == "infinity": 
return ((x_0,0),r) 

else: 
L = P[0O]-r 


m =(0-y_1)/(L-x_1) 

#f (x) =m(x-x_1)+y_1 

#g (x) =x_0 

y_2 = m*(x_0-x_1)+y_1 
return(x_0,y_2) 


We now have the transformed segment AD, and the next step is to find the 


perpendicular bisector of this vertical segment. 


Proposition 5.4. If there exist a circle e with center F = (xo0,0), and radius r, then 


r?=FD-FA. 


Proof. Let A = (x0, yo), D = (#1, y1), and F' = (0,0) using the Definition 3.11 we create 
a circle centered at the origin where the length of the radius takes point D, which in 


Figure 5.10 is inside of circle e to the point A, which is outside of circle e. Similarly for 


A being the inversion of D. So r = y/|yo - y1 — Yo! 


Figure 5.10: Circle e 
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In Sage, I create a function for the circle inversion that sends FD to F'A by 


creating circle e. 


def perp_bisectors_segment_h(P_0,P_1): 


xO = P_O[O] 
yO = P_O[1] 
x1 = P_1[0] 
yi = P_i[1] 


radius_F = sqrt(abs(y0*y1-y0)) 
return((0,0) ,radius_F) 


Proposition 5.5. Inversion of circle e over circle d, takes e to circle e' 


Proof. Let circle d have the center (x0, yO) with radius r and circle e has the left endpoint 
(21, y1) and the right endpoint (x2, y2). By Definition 3.11, we know the left endpoint 
goes to the transformed right endpoint, TR (See Theorem 3.14). If rg < 21, then TR = 
(aa) + xo, and if rp > 71, then TR = — (aan) + 2x9. The right endpoint goes to 


the transformed left endpoint, TL. If x9 < xo, then TL = ( ) +20, and if ro > ro, 


2 
|vo—29| 
then TL = — (ma) + x29. Now once you have the transformed left and right points, 
using Proposition 5.2 we find the geodesic between two points, this gives us e’ (See Figure 


5.11). 


Figure 5.11: Perpendicular bisector of segment AB 


Using Proposition 5.5, I create a Sage function that finds the endpoints of circle 
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def perpbisector_e_transformed(P0O,P1,P2,r): 
#PO is the center of the circle d 


x0 = PO[O] 
yO = PO[1] 
#P1 is the left endpoint of circle e 
x1 = P1[0] 
yi = P1i[1] 
#P2 is the right endpoint of circle e 
x2 = P2[0] 
y2 = P2[1] 
#TL is transformed left point from inversion 2 
if x0<x1: 
TR = (r72/(abs(x0-x1)))+x0 
else: 


TR = -(r°2/ (abs (x0-x1)))+x0 
#TR is transformed right point from inversion 2 
if x0<x2: 

TL 
else: 

TL = -(r°2/ (abs (x0-x2) ))+x0 
return ((TL,O),(TR,O)) 


(r*2/ (abs (x0-x2) ) )+x0 


The created function above requires lots of information: circle d’s center and radius, 
and circle e’s endpoints. I loop it all together to keep track of the orbit point with its 
associated perpendicular bisector, this list is called data3. Recall from Proposition 5.2, 
that if A,B have the same x-values, but different y-values, then the geodesic that goes 
through both points is centered at (29,0) and has a radius of infinity. If this is the case, 
then I do not need to move point b to the vertical line that goes through AA’. I use 
Proposition 5.4, to find circle e. If points A and B have the same y-values, then the code 


uses Proposition 5.4 to find the perpendicular bisector of the vertical segment. 


data3 = [] 
for w in range(1, len(complex_to_real_even_i_points_list)): 
try: 


original_point = (complex_to_real_even_i_points_list[w] [0], 
complex_ to_real_even_i_points_list[w] [1]) 
#pt from orbit 
if original_point[0]== 0: 
perp_bisectors = perp_bisectors_segment_h(original_point, 
(0,2)) 
data3.append([original_point,perp_bisectors] ) 
else: 
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transformed_point = circle_inversion1((0,2), original_point) 
circle_C = hyperbolic_geodesic_circle((0,2), original_point) 
#starting circle 
circle_C_endpoints=circle_endpoints(circle_C[0], circle_C[1]) 
#starting circle endpoints 
left_endpoint_for_circle_C = circle_C_endpoints [0] 
circle_D = radius_of_circle_d(left_endpoint_for_circle_C, 
(0,2)) 
#circle D 
circle_D_radius = circle_D[1] 
circle_D_endpoints=circle_endpoints(circle_D[0], circle_D[1]) 
segment_H = perp_bisectors_segment_h((0,2), 
transformed_point) 
#perpendicular bisectors of a and b’ 
circle_e_endpoints= 
circle_endpoints(segment_H[0] ,segment_H[1]) 
circle_eprime_endpoints= 
perpbisector_e_transformed(circle_D[0], 
circle_e_endpoints[0], 
circle_e_endpoints[1], circle_D_radius) 
perp_bisectors = hyperbolic_geodesic_circle 
(circle_eprime_endpoints [0] ,circle_eprime_endpoints[1]) 
data3.append([original_point, perp_bisectors] ) 

except: 

print ("line") 


This list produces all of the perpendicular bisectors of segments AB. We now 
have enough information to build the Dirichlet domain. I start by graphing the perpen- 


dicular bisectors. 
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Figure 5.12: Graph of all perpendicular bisectors between segment AB 
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Shown in Figure 5.12, there are more walls than necessary, so I clean up and 


remove most of them to get Figure 5.14. 
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Figure 5.13: Hyperbolic polygon 
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5.4 Analzying the Polygon 


I start analyzing this surface by finding the area. The Sage code, data3 produces 
the center and radius of the perpendicular bisectors. Using this information, I find the 


angle of intersection between each circle using Darboux’s Theorem. 


Theorem 5.6. [Wei99, Circle-circle intersection] The angle of intersection between two 


> ae} 
circles is given by cos(@) = a, where d is the distance from the center of the two 


circles, and r1,1r2 are the radii of the two circles. 


This computation produces the external angle of intersection between the two 
circles. There are two intersections where it is the internal angles I actually need to find 


the area. In order to get this internal angle, I compute 7 — cos@. Using Theorem 5.6, I 
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produce angles with radian measures of: 
e 1.69612415796296, 


e 1.69612415796296, 


1.69515132134168, 


2.09439510239319, 


e 1.44644133224814, 


1.44546849562683, 


1.44546849562680, 


e 1.44644133224809, 


2.09439510239319, 
e 1.69515132134167. 


Using the Gauss-Bonnet Theorem, Theorem 2.25, I compute an area of 8.37758040957283 
radians. Using the Arithmetic Area Formula 4.14, I compute an area of “ & 8.3775804096. 
Using the two different area formulas is a way to check that this hyperbolic polygon is 


the 10-gon that I say it is, with the correct walls (i.e. I am not missing any). 


Proposition 5.7. If 7 is a mapping that moves a Dirichlet domain, F', to an adjacent 
domain, then F to an adjacent domain Dirichlet domain with walls W = {w1,..., wr}, 


then yw; € W. 


Proof. If F is a Dirichlet domain with walls W = {wi,...,w,} and y moves F' to an 


adjacent domain, then 7F' has walls {yw,...,yw,}. If FO y7F is w;, then there exist 


j € {1,...,r} such that w; = yw ;. Therefore, y~tw; € W. 


Example 5.8. Let’s look at the complex point (—2.8922386553849+0.442678727629470:), 
which is a point I produce from the orbit sample. The perpendicular bisector that correlates 
to this point is shown in Figure 5.14 below, this is the pink geodesic that has the center 


(—3.71437635490060,0). The generator that corresponds to this point is 1/2 —i—1/2k. 
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The inverse of this generator is 1/2 +%i+1/2k. Looking a Figure 5.14, the perpen- 
dicular bisector that has this generator is also in pink to the left, and has the center 
(—1.99554730922155,0). The orbit point that corresponds to this perpendicular bisector 
is (—1.74554297357772 + 0.2505621735837027). Similary, if I start with the orbit point 
(—1.74554297357772 + 0.2505621735837027) the inverse of the generator corresponds to 
the point (—2.8922386553849 + 0.4426787276294707). It actually does this for every wall 
in the polygon. Using the symmetry of generator mappings, I produce a list and image of 


which walls were glued together, see Figure 5.14 below. 


Figure 5.14: Hyperbolic polygon with gluings 


Notice that the pink walls and the orange walls, are the only walls that are 
adjacent. This is due to them corresponding to elliptic elements, which means that 
Mobius transformation is a rotation (see Proposition 3.10). To confirm these are elliptic 
elements, I computed the trace, and it is 1 for both the pink and orange mappings. A 
rotation produces cone points in the quotient. All of the other mappings are hyperbolic 
elements, which means they have a trace that is greater than 2. These transformations 


are translations. Using these edge gluings, I create a hyperbolic surface. 
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Figure 5.15: Hyperbolic Surface 
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In Figure 5.15, we can see that topologically I produce a genus 1 surface with 
two cone points, i.e. it has signature (1;3,3). To see what the cone points are, I use 


the Riemann-Hurwitz Theorem to verify the area (See Theorem 2.29). By plugging in 0 


punctures, 1 genus, and two cone points of 3,3 we have 2m(2(1) — 2+ 2(45+)) = &. This 


confirms that this hyperbolic surface has two cone points of order 3. 


Theorem 5.9. The hyperbolic orbisurface created from the quaternion algebra (3,5) over 
Q has signature (1;3,3;0), and can be cut into a hyperbolic 10-gon, depicted in Figure 
5.15 and Figure 5.14. Furthermore, the following order elements both glue together the 
walls of the polygon and generate the group OW): 
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1. 1/2 +i + 1/2k, 
2 1/2 42 179k: 
3. -1/2 + i - 1/2k, 
4. 1/2 +4 - 1/2, 
o. 2+ 4, 
6 B24, 
7. 8/2-5 + 1/9, 
8. 38/2 +5 - 1/9, 
9. 3/2 - 5 - 1/9k, 


10. 8/2 +5 + 1/2k 
Corollary 5.10. With the notation as above, if T is the image of O in PSL2(Q{d])), 


d = V3, under the standard embedding (Definition 4.6), then T is generated by the ma- 


trices: 
; s+d sd 
1g lig 
1_q 3d 
9 |2 2 
sd $+ i 


10. 


rr | 

rr | 

rs | 
rs | 
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Chapter 6 


Conclusion 


The polygon I produce is created by performing a careful analysis of the principal 
arithmetic surface (3,5) over Q. To find this hyperbolic polygon, I use Sage to adjoin the 
two symbols 7? = 3, 7? = 5, and produce the quaternion algebra over Q. The quaternion 
algebra produces a list of matricies that I use to transform the point (0,2), and find the 
orbit sampling. Taking each point from the orbit sampling, I create segments between 
them and (0,2). I then find the Dirichlet domain by finding the perpendicular bisectors 
of the segments and look at finitely many, but enough walls. Once I have the Dirichlet 
domain, I carefully analyze this surface by finding the hyperbolic area, the associated 
gluing maps, and the matrices that correspond to gluing one side to the other. 

Creating a hyperbolic polygon is more than just an example problem. From the 
quaternion algebra, we produce a list of generators, and these generators are used to find 
the matrices used to transform the point (0,2). When you find the Dirichlet domain, 
you have a specific list of generators that correspond to side pairings of the walls of the 
hyperbolic polygon. Theorem 5.9 is valuable to other researchers for many reasons. This 
list of finite generators has uses in other areas of math. Researchers use the generators in 
systolic geometry, which is the study of invariants of manifolds and polyhedra to find the 
systole of the surface and its covers. It can be used to find the Cheeger constant of the 
surface, which is a numerical measure of whether or not the surface has a “bottle neck”. 
This data is useful for checking future results and hypotheses. 

The process leading up to finding my hyperbolic polygon leaves me with a few 


unanswered questions. It is an open question as to how Macasieb excludes certain sig- 
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natures from her list of genus 2 surface coverings. I am unable to construct a few of 
the signatures on her list. It is an open question as to how they can be constructed or 
if they cannot and should be excluded from the list as well. Generalizing this to the 
genus 3 surface, is it possible to construct the other signatures Proposition 2.34 for a 
genus 3 surface covering. Once I have these coverings how do I confirm that this list is 
correct. Next I want to reevaluate {[Mac08] and surface coverings. I really enjoy making 
the constructions of these surface coverings because I feel like it pushes me outside of my 
comfort zone. I have never considered myself as a creative person, but in order to make 
these constructions, I have to be creative in the bending and stretching of these surfaces, 
as well as where to place the axis of rotation. I want to continue to challenge myself in 
this area and advance myself as a student and person. 

This thesis brought together many areas of math, including topology, hyperbolic 
geometry, linear algebra, group theory, and noncommutative algebra. Additionally this 
project involves a significant programming component, in which I write and implement 
code in the Sage. My project has expanded my mathematical thinking by challenging me 
to expand my understanding of the algebra that is behind the Euclidean and hyperbolic 
geometry of these surfaces and figures, like GeoGebra images or surface constructions, 
and put it into a Python-based Sage code. With surfaces, I have to challenge my mind to 
see the bending, stretching, and gluings of surfaces. For example, taking my hyperbolic 
polygon, and finding its corresponding surface (see Figure 5.15). I have never consider 
myself as a creative person or visual learner, but I feel this thesis challenges me the most 
in this aspect and make me a better learner. Going into this thesis, I have minimal 
experience with Sage and I have never created my own functions. I found it extremely 
difficult to think like a “coder”. I often use GeoGebra to create the images for which I 
need to create functions and piece together how to get Sage to output the image I am 
seeing. This is a process. By the end of my thesis, I am able to create my own functions, 
pull out the data that I need from certain list, and create images (see Figure 5.14). This 
thesis has expanded my learning by helping me to step outside of my comfort zone, and 
change my way of thinking as well as challenging me to bring together multiple topics 
such as topology, hyperbolic geometry, linear algebra, group theory, and noncommutative 


algebra. 
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